Analytical approach to model of scientific revolutions 
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The model of scientific paradigms spreading throughout the community of agents with memory 
is analyzed using the master equation. The case of two competing ideas is considered for various 
networks of interactions, including agents placed at Erdos-Renyi graphs or complete graphs. The 
pace of adopting a new idea by a community is analyzed, along with the distribution of periods 
after which a new idea replaces the old one. The approach is extended for the chain topology onto 
the more general case when more than two ideas compete. Our analytical results are in agreement 
with numerical simulations. 
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I. INTRODUCTION 

There is a tendency to separate certain periods in the 
history of civilizations, such as Renaissance or Enlight- 
enment, which qualitatively differ from each other by 
dominating trends in science, art or customs. Techno- 
logical innovations and scientific discoveries constantly 
emerge and it is not likely that some kind of equilib- 
rium — "end of history" [l| — will ever be reached. 
The changes (evolutionary and revolutionary ones) hap- 
pen due to the interactions and exchange of innovative 
ideas [2.-^] at the level of individuals, communities or 
even civilizations. Eventually, ideas spread throughout 
the communities d, 01- Some of the ideas gain broad 
(even global) acceptance and popularity, replacing old 
ones Similar phenomena can be observed in models 
of opinion formation dynamics l9l-fl2j . 

The process of adoption of an innovative technology 
or a new scientific concept by individuals and com- 
munities differs from adoption of, for example, a new 
trend in arts. Obsolete technologies and discarded sci- 
entific theories, once abandoned, are not likely to be ac- 
cepted by individuals again. To model such a process, 
agents should be given some kind of memory. Another 
important fact is that the will of individuals to adopt a 
new scientific concept depends on its global popularity. 
For example, spreading of technological innovations is 
often slowed down by incompatibility with existing stan- 
dards. In the field of arts, the situation is different. Old 
ideas can reemerge and become popular again, such as 
Renaissance artists were inspired by Antique philosophy 
or architecture. 

Recently a model was introduced by Bornholdt et al. 
which attempts to describe scientific revolutions (TJ). 
The model combines interactions at the level of individ- 
uals with infiuence of the whole community. Despite its 
simplicity, it managed to reconstruct some key features of 
the dynamics of scientific paradigms spreading, including 
an asymmetry between the rate of adopting a new idea 
by the community and the speed of its decline when new 
competing ideas appear. The model was based on nu- 
merical simulations and no analytical treatment was pre- 



sented. In this paper, the master equation and Markov 
processes theory [15| were applied to analyze the dynam- 
ics of the system in the case of small level of agents' cre- 
ativity for various topologies of agent interactions, such 
as chain, complete graph, ER graphs, star and square 
lattice. For chain topology the approach was extended 
in an attempt to describe the system dynamics for higher 
levels of creativity. 



II. THE MODEL OF SPREADING OF IDEAS 

The rules of the model [ij are very simple. N agents 
occupy nodes of a network. Every agent follows some 
paradigm (idea), labeled by a natural number. In each 
time step a random agent i (with paradigm Si) is se- 
lected, along with one of its neighbors j (with paradigm 
Sj). If the agent i has never followed the paradigm Sj, 
it adopts it with probability NgjN, where Ng- denotes 
the number of agents representing paradigm Sj. Addi- 
tionally, new paradigms, which have never been present 
in the community, can appear: with probability a a ran- 
dom agent is selected, which changes its paradigm into 
one which has never been present in the community. 

The most important feature of the model is the mem- 
ory of the agents, who do not adopt the same paradigm 
twice. One can find analogy between this model and evo- 
lutionary dynamics models: innovations can be regarded 
as mutations which allow the affected individuals outper- 
form their rivals. Lack of any evident fitness parameter, 
which would describe how well a specie is adapted to the 
environment, is not necessarily a drawback of such an in- 
terpretation, as the fitness of a specie is alway a posteriori 
knowledge [16]. 

There are some general features of the evolution of the 
system, independent of the interactions network topol- 
ogy. For very small probability a at most two paradigms 
coexist (other cases are neglectable due to their much 
smaller probability). This case will be analyzed for vari- 
ous networks in Sec. IIIIl 

For higher values of probability a, other effects have 
to be considered. In such a case usually more than two 
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paradigms coexist, which "compete" with each other. 
However, one may suppose that within a relatively wide 
range of a still two paradigms can be separated at every 
moment: the "old" paradigm which is the most popular, 
but currently at the decline, and the "new" paradigm, 
the second most popular, which will prevail after some 
time (and then enter the stage of decline). This case will 
be analyzed in Sec. IIVI 



III. THE CASE OF TWO COMPETING 
PARADIGMS 



new paradigm), ki is the degree of node i and Uij is the 
adjacency matrix. For n = N so defined transition rate 
is automatically equal to 0, since V^Si = then. 

It can be easily proved that if all the transition rates 
are different (k ^ j ^ Wk ^ Wj), the solution of Eq. ^ 
with the initial condition P{n, 0) — dni is 



P(n,t) = ^Cre-^'=*, 

k=l 



(5) 



where 



A. General case 

When the creativity level of the agents a is small 
enough, at most 2 paradigms coexist, referred to as the 
paradigm (at the stage of decline) and the paradigm 
1 (at the stage of expansion). The evolution consists of 
two distinct periods. 

1. All the agents share the same paradigm 0. The 
length of this stage of stagnation is a random vari- 
able of the exponential distribution 



(1) 



and the mean value 



stag / 



Tstag=0 

1 - a 1 



a a 



(2) 



2. After an innovative paradigm 1 appears, it starts 
spreading across the community. The time of ex- 
pansion of the paradigm 1 will be denoted T. It is 
a random variable whose distribution depends on 
the interactions network topology. After time T all 
the agents share the paradigm 1 and the state of 
the system is equivalent to the initial one. 

In our approach the state of the system is character- 
ized by one variable — the number n of agents shar- 
ing paradigm 1. The problem reduces to the problem 
of expansion of the paradigm 1 throughout the commu- 
nity, starting from one agent with the paradigm 1 at time 
t = 0. The generic master equation has only two terms: 



^P(n,t) = P(n-l,t)VF„,„_i 



P{n,t)Wr 



n+l,n 7 



(3) 



where transition rates from the state n to n + 1 (for n G 
[1, — 1]) are equal to 



n 

iV2 



^ 1 



i=l 



k. 



N 

-E 



(4) 



where Si denotes the state of ith agent (s^ = ^ agent 
follows the old paradigm, = 1 agent follows the 



n— i n 

c?. = T\WiT\ . 



(6) 



Let us note that since — and Vi<„<ArW„ > 0, 
the distribution evolves into 



lim P(n, t) = 6n 



N 



(7) 



(all the agents share the paradigm 1), which is an ex- 
pected limit. 

The approximation of only two competing paradigms 
makes sense if the mean stagnation time (Tstag) = l/a is 
greater than the mean expansion time (T), i.e. 



a < 



1 



(8) 



This upper limit of a has to be estimated for each type 
of network separately. In general the P{T) distribution 
can be expressed by the P{n, t) probability as 

P(T = t) = WN-iPiN - 1, t - 1) « ^PiN - 1, t) (9) 
and, taking into account Eq. ([5]), we obtain 



1 %=^' -1 



tPiT^t)dt.-Y^^^^ 



(10) 



B. Chain topology 

Let us consider the case when the agents occupy nodes 
of a chain. For simplicity, periodic boundary conditions 
will be assumed. 

This specific topology makes the problem quite sim- 
ple and in the first approximation it is not necessary to 
analyze the master equation (|3]). The average number 
of agents sharing paradigm 1 can be derived from the 
recursive equation 



(n(0)) = 1 
(n(t + l)) = 



(n(t)) 



2 1 ("(«)) 
N 2 N 



l{t)) (1 



1 



(11) 
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which has the sokition 



On the average, after time 

(T) = log N 



(12) 



(13) 



the paradigm 1 will stop spreading as it will be shared by 
the whole community. The situation will be stable un- 
til another innovation appears. From this condition the 
range of a can be estimated for which this approximation 
makes sense. Taking into account Eq. ([8]) we get 



a < 



1 



7V2 logiV 



(14) 



For a more exact analysis one has to consider the 
master equation Let us make a simple observation 
that the subgraph consisting of agents sharing the new 
paradigm is connected. The transition rates are therefore 
equal to 



^ 1 

1=1 



(15) 



Above we used the periodic boundary conditions, so 
agents at positions 1 and N + 1 are equivalent. 

To solve the problem, at first we will neglect the 5nN 
term and treat the n variable as if it could grow to in- 
finity, n = 1,2,. ..,c». Eventually, the transition rates 
from state n to n -I- 1 are equal to Wn = jp- and the 
master equation has the form 

-P(n,i) = P{n-l,t)^ PM^- (16) 

Due to the simple form of the transition rates, Eq. 
can be solved using the method of characteristic function 
G: 



G(s,t)EE(e™«). 



(17) 



This approach has such an advantage over using ([S|), that 
solutions are automatically in a compact form. The mas- 
ter equation (|T6)) with the initial condition P{n, 0) = Sno 
leads to the partial differential equation with the initial 
condition 

fG(.,t) + ^(e--l)^G(s,0=0 
G(s,0)=e*^ 



which can be solved as 

G{s,t) = 



1 



After a short algebra it can be proved that 



(19) 
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FIG. 1: (Color online) Chain graph topology, A'^ = 64 nodes, 
a < Evolution of the system starting from P(n,t = 

0) — Sni; probability P{n,t) at various moments t. Points are 
obtained from the numerical solution of the master equation 
(fT6|) . Lines — analytical predictions (Eq. (|23|) \ 



so 



P(n,t) = e-*/^' (l-e-*/^') 



n-l 



(21) 



This is valid for n < A^. In order to take into consid- 
eration the limitation on the n variable (the SnN term 
in Eq. ([T5|)). one has to consider the accumulation of 
probability at point n = N: 



m — 1 



1 - e~*/^' 



N-l 



(22) 



Eventually, 



(23) 

The mean value of n resulting from the distribution ([25)1 
is equal to 



n) = e*/^' ( 1 - ( 1 - e 



-t/N- 



N 



(24) 



which for small t reduces to ([1] 

From the distribution P{n,t), a more exact approxi- 
mation of (T) than Eq. can be obtained. According 
toEq. ®, 



ri,=l 



P(T = i)«le-/^^(l-e-/^y"', (25) 
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FIG. 2: (Color online) Chain graph topology, = 64, 128, 256 
nodes, a < 1/{T). Evolution of the system starting from 
n — 1 innovative agent: in) versus time. Points — simulated 
data. Lines — analytical predictions (Eq. (|24|) '). 
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FIG. 3: (Color online) Chain graph topology, a < 1/{T). 
Distribution of the expansion periods lengths T for different 
system sizes TV. Points — simulations, lines — analytical 
predictions (Eq. (|25|) 'l. (Inset) Mean time of expansion (T) 
versus system size TV; simulations (points with error bars) 
compared with the analytical predictions (line) (Eq. (|26p ). 



(T) = 



Y^tP{T = t) 

t=Q 



te 



JV-2 



dt 



7V-1, 



(26) 



where is the nth. harmonic number. Since harmonic 
numbers grow approximately as fast as the natural log- 
arithm, the approximated solution (I13|) is very close to 
dill). In fact, for TV > I 



Hn ~ logTV + 7, 



(27) 



where 7 « 0.5772 denotes the Euler-Mascheroni con- 
stant. 



C. Complete graph topology 

Let us consider the situation when interaction is pos- 
sible between every pair of agents, i.e. V(i j)aij = 1. Re- 
ferring to the generic master equation (|3]), the transition 
rates are equal to 



N N 

n \ - 1 ~ Si X - 

TV2 ^ TV - I '"^ 



TV2(TV- I) 



TV3 



Eventually, the master equation has the form 

d , ,(n- l)2(TV-n + l) 



(28) 



TV3 



n^{N - n) 
TV3 ■ 



(29) 



If all the transition rates are different {j ^ k ^ Wj 7^ 
Wk, which is satisfied if equation N = a{l + b + b'^) does 
not have trivial solutions a, 6 among natural numbers), 
the solution can be written in the form of the sum ([S|): 



(30) 



fc=i 



where 



3^k 



\P{N-j)~kHN-k) 



1 



(TV - n)! p{N - j) - k^{N 



The expansion time T distribution P{T) can be derived 
from the P{n, t) distribution: 



N-l 



P(T ^t)^ ^P{N~l, t) = ^ E ^k-'e 



k=l 



(32) 



The analytical predictions are in agreement with the 
simulations (Fig. HHSl)- 



D. Erdos-Renyi graph topology 

Let us consider the situation when the network of in- 
teractions is an Erdos-Renyi graph . For each pair of 
nodes the edge between them exists with probability p. 
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FIG. 4: (Color online) Complete graph topology, A'^ = 64 
nodes, a < 1/{T). Evolution of the system starting from 
P{n,t = 0) = Sni; probability P{n,t) at various moments t. 
Points are obtained from the numerical solution of the master 
equation ([29|) . Lines — analytical predictions (Eq. (f30|) l. 



FIG. 6: (Color online) Complete graph topology, a < 1/{T). 
Distribution of the expansion period lengths T for different 
system sizes TV. Points — simulations, lines — analytical 
predictions (Eq. ((32)) '). 
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FIG. 5: (Color online) Complete graph topology, — 
32,64, 128 nodes, a < 1/{T). Evolution of the system start- 
ing from n = 1 innovative agent: (n) versus time. Points 
— simulated data. Lines — analytical predictions (obtained 
from Eq. (O). 



The degree distribution of an ER graph is a binomial 
distribution 



P(fc) = Bin{N,p) 



N-1 



(33) 



The transition rates of the master equation ^ are 
equal to 



i=l 3 = 1 



N 



N 



N 



7V2 



i=l 



kj + k; 



(34) 

where the random variable denotes number of i's 
neighbors following the paradigm 1 and k~ — the 



paradigm 0. Within the mean field approach, the transi- 
tion rates can be estimated by 



n{N — n) 



fc+ 



k+ + k- 



(35) 



where (•) denotes averaging on the whole population of 
the agents. Treating A:"*" and k~ as independent random 
variables with binomial distributions, we obtain 



k+ 



k+ + k- 



E E 

k+ = l k-=Q 



E^^ 



k+ 



k+ + t 



-p{k+)p{k~) 



fc+=i 

\ k- y„fe"/i „\N—n-l-k- 



E 



k+ + k 



^p" {1-pY 



k+ = l 



k+ 



l_p\ J ^fc+-i(i+^)^— irf^ 



N- 



-(l-(l-p)^-i) 



(36) 



Eventually, 



n'^{N - n) 



(l_(l_p)A^-l) 



(37) 



For p = 1 the transition rates reduce, as expected, to the 
ones obtained for the complete graph topology. 
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FIG. 7: (Color online) ER graph topology, TV = 128 nodes, 
p = 0.5, 0.05, 0.02, Q < 1/(T). Evolution of the system start- 
ing from n = 1 innovative agent: (n) versus time. Points 
— simulated data. Lines — analytical predictions (obtained 
fromEq. 



The master equation has the form 



- P{n,t) 



n'^{N - n) 



[i-{l-p) 



N-l\ 



Similarly to the case of complete graph topology, the so- 
lution can be written in the form of the sum ([5]): 



where 



k=l 



(n 



{N - ny. -^A p{N - j) - k^{N - 



Now the expansion time distribution P{T) is 
P{T = t) « lp(7V_i,i) 

^ _L ^ C'A'-lg-fc'(^V-fc)(l-(l-p)"-')*/^^|l) 
^ k=l 

As can be seen in Fig. [7l our approach predicts a 
decline in the rate of growth of the new paradigm cluster 
with decreasing the network density (parameter p), but 
seriously underestimates that decline. We suppose that 
there are some nontrivial correlations between the agents' 
states resulting from the dynamics, which were not taken 
into consideration. 



E. Star topology 

1. Central agent innovative 

Let us consider the star topology of interactions, i.e. 
there exists a central agent connected to all the other N — 
1 periferal agents — and those are the only connections. 
Moreover, we will require that at time t = the central 
agent follows the innovative paradigm 1. This case is 
interesting since in this variant the new idea spreads at 
the highest pace — the transition rates 

{N - n)n 



TV2 



(42) 



are the highest possible for this model among all the pos- 
sible topologies (one can for example compare Eq. 

to Eq. (HSJ, (1211), (ESD). 

The periferal agents are not connected with each other 
and they are only influenced by the "mean field" of 
paradigms. Therefore, one can consider the change in 
time of the average state of a periferal agent. Let Tr{t) 
denote probability that at time t a periferal agent fol- 
lows the paradigm 1. The evolution of 7r(t) follows the 
recursive equation 

7r(0) = 
nit + l) = 4t) + (l-4i))^inM 

= 7r(t) + i(l-7r(i))(^(t) + ^ 



which can be solved in the approximation of continuous 
time: 



1 



(7V-l)exp(-^) + l N-\ 



(44) 
(45) 



Thus, the average number of agents sharing the paradigm 
1 is equal to 



{n{t)) = {N-l)^{t) + l 



N 



1 + iVe-*/JV • 



(46) 



From 7r(t), the expansion time distribution P{T) can 
be derived: 



TVe*/^(e*/^ - 1)^-2 



N 



(47) 



We were not able to find an analytical formula for (T) . 
However, since P(T) is unimodal with a well defined max- 
imum, one can assume that (T) grows as fast with N as 



argmaxP(r) = A^log 



N^-N 



0{N log N) 



(48) 
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FIG. 8: (Color online) Star lattice topology with central agent 
innovative, iV = 64,128,256,512,1024 nodes, a < 1/{T). 
Evolution of the system starting from n = 1 innovative agent 
(central): (n) versus time. Points — simulated data. Lines 
— analytical predictions (Eq. (|46|) 1 



FIG. 10: (Color online) Star topology with all agents equally 
innovative, a < 1/{T). Distribution of the expansion period 
lengths T for different system sizes A^. Points — simulations, 
lines — analytical predictions (Eq. ((51}). 
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FIG. 9: (Color online) Star topology with central agent in- 
novative, a < 1/{T). Distribution of the expansion period 
lengths T for different system sizes A'^. Points — simulations, 
lines — analytical predictions (Eq. (147^ 1. 



Indeed, as can be seen in Fig. the mean expansion 
time can be very well approximated by 



(T) w kN log N, 



(49) 



The system is similar to the one described above. The 
only difference is that at time t = any agent can be 
the innovative one. With probability (A^ — l)/iV 1 a 
periferal agent will become the innovative one and after 
time To, which is a random variable with exponential 
probability distribution 



P{To^t) = 1- 



1 



Ar3 



1 



I) J N^{N-l) 



(50) 



the central agent will adopt idea 1. Then the dynam- 
ics will be as described in the case of innovative central 
agent. Thus, the probability distribution of the expan- 
sion time can be well approximated by a convolution of 
two probability distributions: 



P(T = t) 





1 



N 



N3 



t). 



(51) 



The mean and the standard deviation of this exponen- 
tial probability distribution are equal to 



(52) 



As can be seen in Fig. [TTl the analytical results for 
both variants of star topology agree with simulations very 
well. 



where parameter k, as obtained by fitting to simulated 
data, is equal to A: = 2.149 ± 0.007. Since the transition 
rates (I42p are the highest possible, the mean expansion 
time must be the shortest for this model among all 
the possible topologies. 



F. Square lattice topology 

Let us consider the square lattice topology. Periodic 
boundary conditions will be assumed, so each agent has 
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FIG. 11; (Color online) Star lattice topology. Mean time 
of expansion (T) versus system size A''. Empty circles — 
innovative central agent, filled circles — all agents equally 
innovative. Error bars correspond to the standard deviations 
of the samples. Lines — analytical predictions, Eq. (|49[) and 
(EH). 



FIG. 12: (Color online) Square lattice topology, = 
64, 169, 256 nodes, a < 1/{T). Evolution of the system start- 
ing from n = 1 innovative agent: (n) versus time. Points 
— simulated data. Lines — analytical predictions (obtained 
from Eq. ([5]) with transition rates ((53])). 



4 neighbors. The first approximation would be to assume 
that the cluster of agents sharing paradigm 1 grows uni- 
formly in each direction, so at any moment it is circle- 
shaped, with the radius of the circle equal to r = ^njix. 
Therefore, within this approximation the transition rates 
in the generic master equation ([3]) are equal to 



,3/2 



-(1 - <5„Ar). 



(53) 



Similarly as in the case of complete graph topology, 
the solution of the master equation 



d 



fPin^t) = P{n~\,t) ^^''^J^ (l-<5„,i,Ar) 



(54) 



can be expressed in the form of the sum of products ([S]) . 
Having compared the results of such an approximation 
with the simulations (Fig. [T^ . one has to state that 
this approach significantly overestimates the pace of the 
growth of the new paradigm cluster. 



IV. THE CASE OF MANY COMPETING IDEAS 



our results for chain topology, which is probably the sim- 
plest one. 

For a higher than 1 / iV^ log N Eq. ([TT|) has to be ex- 
tended by terms describing the appearing of new clusters 
of ideas, which would slower the process of expansion of 
the paradigm 1. In the first approximation it will be 
assumed that the only important new clusters are those 
appearing inside the cluster of the paradigm 1 and that 
they do not overlap with each other. Their growth is 
described by Eq. (fT2|). Thus, the recursive equation for 
{n{t)) is now 



{nit + 1)) = (nit)) ( 1 + 



1 



^ (n(T)) 

2^ « ^ (An„eto (t~T)} 



-^^(n(T))exp 



(56) 



Substituting the sum with the integral and stating that 
{n{t)) — exp {t/N'^) f{t) leads to the following equation 
for fit): 



If the mean stagnation time {Tgtag) 
than the mean expansion time (T), i.e. 



a = 



> 



stag / 



(T)' 



1/a is shorter 



(55) 



it is most probable that more than two paradigms coex- 
ist in the community at any moment. This case is much 
more difficult to describe analytically. Below we present 



a 



(57) 



which, assuming the same initial conditions as in Eq. 
([Tl]) (single innovation at time t = 0), has the solution 



fit)=cosiXt), 



(58) 



where A = y^a/N^. The complete formula for the first 



9 



approximation of {n{t)) is therefore 



{n{t)) exp ( — ) cos (Xt) . 



(59) 



As it should have been expected, for a — > this approx- 
imation converges to the previous one (Eq. ([T2|)). 

A better approximation can be obtained by substitut- 
ing the term exp (^^) by {n{t — r)) in Eq. ([56]), as new 
paradigms can also be "attacked" by paradigms appear- 
ing after them. The equation 



(60) 



(n(i + l)) = {n{t)) [l + j^ 



T = 

does not have a simple analytical solution, but by sub- 
stituting the sum with the integral and stating {'n{t)) = 
exp (t/N^) cos(At)(l -I- g{t)), where g{t) <^ 1, an integral 
equation can be obtained, 



= -Asin(A<) +.g'(t)cos(Ai) 



(61) 



+ \' / cos(At)(1 +.g(r)) 

• cos(A(f - t))(1 + g{t - T))dT, 

which can be is solved provided all the terms in the inte- 
gral apart from the product cos(Ar) cos(A(i — r)) are ne- 
glected. Eventually, the second approximation of (f^(^)) 
obtains the form 

{n{t)) = exp j cos(At)(l - log | cos(Ai)| - hh^). 

(62) 

The comparison with the simulations (Fig. I13p shows 
that the last approximation (Eq. (|62|) ) is better than the 
previous ones (Eq. (HI]), Eq. ((5^ ). 



V. CONCLUSIONS 

We have developed an analytical approach based on 
master equation that describes a model of paradigms evo- 
lution |14] and confronted our results with the outcome of 
the work of Bornholdt et al. as well as with our numerical 
simulations. The outcome suggests, that the asymmetry 
between the paces of growth and decline of the dominant 
idea, observed in [3] , is a generic property of the model 
and should be observed for any topology of interactions. 

Our analytical methodology can be used to consider 
various topologies of interaction networks. The crucial 
parameter of the dynamics is the creativity of the agents, 
described by the a parameter. In the case when agents 
are almost non-innovative, the evolution consists of sub- 
sequent periods of stagnation (single paradigm present 
in the community) and periods of expansion (an innova- 
tive paradigm spreads across the community replacing 



the old one). The mean length of the stagnation pe- 
riod is equal to (Tstag) = regardless of the inter- 
action network topology. On the other hand, the mean 
length of the expansion period (T) strongly depends on 
the topology. If a ^ {T)~-^, the mean time between 
shifts of dominant paradigms can be approximated by 
{Tstag + T) ~ (Tstag) = which is the scaling ob 



served in the simulated data by Bornholdt et al. 
(note that a different time scale was used in [l4|). 

Our approach is mainly based on the approximation 
of two competing paradigms, which is justified if the 
level of creativity a is small enough, i.e. a < 1/{T). 
For each type of interaction network topology this range 
has to be calculated separately. Five different topologies 
were taken into consideration: chain, complete graph, 
ER graphs, star and square lattice. 

For the chain topology it was possible to find compact 
forms of the analytical solutions. The mean expansion 
time (T) scales with the system size N like iV^ log N and 
during the stage of expansion^ the mean size of the cluster 
of the new idea grows like a damped exponential function, 
see Eq. (HH). The analytical results are in agreement 
with the simulated data. 

For the complete graph topology the proposed ap- 
proach also results in a good agreement with the sim- 
ulations, but compact forms of the functions describing 
the system evolution (Eq. (150]) . Eq. ([5^ ) probably do 
not exist. 

In the case of ER graphs topology, our approach under- 
estimates the rate of the growth of the new idea cluster 
(Fig. [7]). The reason can be correlations between the 
node degree and the state of the agent located at that 
node, which were not taken into account. 

Comparison of two variants of star topology brought 
interesting results. In the first variant, when we require 
that the innovation firstly appears in the central node, 
the rate of the expansion of the new idea is the fastest 
possible among all the topologies and the mean expan- 
sion time (T) (which is the shortest possible) scales like 
TV log TV. However, if we remove this requirement and let 
the innovation appear in any node with equal probabil- 
ity, (T) grows to (higher than at a chain, where mean 
distance between nodes is much larger) and almost whole 
that time is taken by "convincing" the central agent to 
the new idea. 

In the case of the square lattice, the method only qual- 
itatively reproduces the results of the simulations (Fig. 
[T2|) . The problem probably lies in the apparently too 
rough estimation of the shape of the cluster of the new 
idea as a circle. 

For a higher level of creativity a, when most of the time 
more than two ideas coexist, the dynamics of the system 
can be found starting from the results obtained for the 
case of lower levels of a and using the method similar 
to the perturbation method. This approach proved to be 
useful in the simplest case — the chain topology. We were 
considering the function describing the mean number of 
agents following the expanding paradigm 1. It was found 
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FIG. 13: (Color online) Chain topology, = 128, various levels of creativity a. Comparison of approximations Eq. (|12p . Eq. 
(|59[) and Eq. (|62p . The red crosses refer to the simulated data. 



that the unperturbed function 

{n{t)) = e*/^^ 
should be modified by two factors. The first one, 



/ a 



< 1, 



(63) 



(64) 



describes "attacking" of the paradigm 1 by paradigms 
appearing after it. The second one, 



1-loe 



1 



-f— V 

4 \N^) 



> 1, 



(65) 

describes "attacking" of the paradigms "attacking" 
paradigm 1. Both these terms, as expected, converge 
to 1 if the creativity of the agents a converges to 0. 

Our analytical approach allows for a better under- 
standing of the system dynamics described by the model 
[l^ and explains some of relationships previously ob- 
served in the simulated data. The proposed methodol- 
ogy can be used to analyze the dynamics of paradigms 



spreading in other networks [T3l • It is especially interest- 
ing since the real networks of human contacts (including 
scientific collaboration networks) exhibit some nontrivial 
properties, such as scale- free behavior [l8|. Investigations 
of such networks are planned in the future. 
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